The evolution of transposable elements in Brachypodium distachyon is governed by purifying selection, while neutral and adaptive processes play a minor role.

Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.

Three near-complete genome assemblies reveal substantial centromere dynamics from diploid to tetraploid in Brachypodium genus.

BACKGROUND: Centromeres are critical for maintaining genomic stability in eukaryotes, and their turnover shapes genome architectures and drives karyotype evolution. However, the co-evolution of centromeres from different species in allopolyploids over millions of years remains largely unknown. RESULTS: Here, we generate three near-complete genome assemblies, a tetraploid Brachypodium hybridum and its two diploid ancestors, Brachypodium distachyon and Brachypodium stacei. We detect high degrees of sequence, structural, and epigenetic variations of centromeres at base-pair resolution between closely related Brachypodium genomes, indicating the appearance and accumulation of species-specific centromere repeats from a common origin during evolution. We also find that centromere homogenization is accompanied by local satellite repeats bursting and retrotransposon purging, and the frequency of retrotransposon invasions drives the degree of interspecies centromere diversification. We further investigate the dynamics of centromeres during alloploidization process, and find that dramatic genetics and epigenetics architecture variations are associated with the turnover of centromeres between homologous chromosomal pairs from diploid to tetraploid. Additionally, our pangenomes analysis reveals the ongoing variations of satellite repeats and stable evolutionary homeostasis within centromeres among individuals of each Brachypodium genome with different polyploidy levels. CONCLUSIONS: Our results provide unprecedented information on the genomic, epigenomic, and functional diversity of highly repetitive DNA between closely related species and their allopolyploid genomes at both coarse and fine scale.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Group VII ethylene response factors forming distinct regulatory loops mediate submergence responses.

Group VII ethylene response factors (ERFVIIs), whose stability is oxygen concentration-dependent, play key roles in regulating hypoxia response genes in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) during submergence. To understand the evolution of flooding tolerance in cereal crops, we evaluated whether Brachypodium distachyon ERFVII genes (BdERFVIIs) are related to submergence tolerance. We found that three BdERFVIIs, BdERF108, BdERF018, and BdERF961, form a feedback regulatory loop to mediate downstream responses. BdERF108 and BdERF018 activated the expression of BdERF961 and PHYTOGLOBIN 1 (PGB1), which promoted nitric oxide turnover and preserved ERFVII protein stability. The activation of PGB1 was subsequently counteracted by increased BdERF961 accumulation through negative feedback regulation. Interestingly, we found that OsERF67, the orthologue of BdERF961 in rice, activated PHYTOGLOBIN (OsHB2) expression and formed distinct regulatory loops during submergence. Overall, the divergent regulatory mechanisms exhibited by orthologs collectively offer perspectives for the development of submergence-tolerant crops.

Systematic analysis of Prx genes in the Brachypodium genus and their expression pattern under abiotic constraints.

Peroxiredoxins (Prx) are ubiquitous peroxidases required for the removal of excess free radicals produced under stress conditions. Peroxiredoxin genes (Prx) in the Brachypodium genus were identified using bioinformatics tools and their expression profiles were determined under abiotic stress using RT-qPCR. The promoter regions of Prx genes contain several cis-acting elements related to stress response. In silico expression analysis showed that B. distachyon Prx genes (BdPrx) are tissue specific. RT-qPCR analysis revealed their differential expression when exposed to salt or PEG-induced dehydration stress. In addition, the upregulation of BdPrx genes was accompanied by accumulation of H2 O2 . Exogenous application of H2 O2 induced expression of almost all BdPrx genes. The identified molecular interaction network indicated that Prx proteins may contribute to abiotic stress tolerance by regulating key enzymes involved in lignin biosynthesis. Overall, our findings suggest the potential role of Prx genes in abiotic stress tolerance and lay the foundation for future functional analyses aiming to engineer genetically improved cereal lines.

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena.

BACKGROUND: The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes). RESULTS: Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa. CONCLUSIONS: The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum.

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

INDETERMINATE1-mediated expression of FT family genes is required for proper timing of flowering in Brachypodium distachyon.

The transition to flowering is a major developmental switch in plants. In many temperate grasses, perception of indicators of seasonal change, such as changing day-length and temperature, leads to expression of FLOWERING LOCUS T1 (FT1) and FT-Like (FTL) genes that are essential for promoting the transition to flowering. However, little is known about the upstream regulators of FT1 and FTL genes in temperate grasses. Here, we characterize the monocot-specific gene INDETERMINATE1 (BdID1) in Brachypodium distachyon and demonstrate that BdID1 is a regulator of FT family genes. Mutations in ID1 impact the ability of the short-day (SD) vernalization, cold vernalization, and long-day (LD) photoperiod pathways to induce certain FTL genes. BdID1 is required for upregulation of FTL9 (FT-LIKE9) expression by the SD vernalization pathway, and overexpression of FTL9 in an id1 background can partially restore the delayed flowering phenotype of id1. We show that BdID1 binds in vitro to the promoter region of FTL genes suggesting that ID1 directly activates FTL expression. Transcriptome analysis shows that BdID1 is required for FT1, FT2, FTL12, and FTL13 expression under inductive LD photoperiods, indicating that BdID1 is a regulator of the FT gene family. Moreover, overexpression of FT1 in the id1 background results in rapid flowering similar to overexpressing FT1 in the wild type, demonstrating that BdID1 is upstream of FT family genes. Interestingly, ID1 negatively regulates a previously uncharacterized FTL gene, FTL4, and we show that FTL4 is a repressor of flowering. Thus, BdID1 is critical for proper timing of flowering in temperate grasses.

Phytochromes transmit photoperiod information via the evening complex in Brachypodium.

BACKGROUND: Daylength is a key seasonal cue for animals and plants. In cereals, photoperiodic responses are a major adaptive trait, and alleles of clock genes such as PHOTOPERIOD1 (PPD1) and EARLY FLOWERING3 (ELF3) have been selected for in adapting barley and wheat to northern latitudes. How monocot plants sense photoperiod and integrate this information into growth and development is not well understood. RESULTS: We find that phytochrome C (PHYC) is essential for flowering in Brachypodium distachyon. Conversely, ELF3 acts as a floral repressor and elf3 mutants display a constitutive long day phenotype and transcriptome. We find that ELF3 and PHYC occur in a common complex. ELF3 associates with the promoters of a number of conserved regulators of flowering, including PPD1 and VRN1. Consistent with observations in barley, we are able to show that PPD1 overexpression accelerates flowering in short days and is necessary for rapid flowering in response to long days. PHYC is in the active Pfr state at the end of the day, but we observe it undergoes dark reversion over the course of the night. CONCLUSIONS: We propose that PHYC acts as a molecular timer and communicates information on night-length to the circadian clock via ELF3.

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei.

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

REGULATOR OF FLOWERING AND STRESS manipulates stomatal density and size in Brachypodium.

Stomata are crucial for gas exchange and water evaporation, and environmental stimuli influence their density (SD) and size (SS). Although genes and mechanisms underlying stomatal development have been elucidated, stress-responsive regulators of SD and SS are less well-known. Previous studies have shown that the stress-inducible Brachypodium RFS (REGULATOR OF FLOWERING AND STRESS, BdRFS) gene affects heading time and enhances drought tolerance by reducing leaf water loss. Here, we report that overexpression lines (OXs) of BdRFS have reduced SD and increased SS, regardless of soil water status. Furthermore, biomass and plant water content of OXs were significantly increased compared to wild type. CRISPR/Cas9-mediated BdRFS knockout mutant (KO) exhibited the opposite stomatal characteristics and biomass changes. Reverse transcription-quantitative polymerase chain reaction analysis revealed that expression of BdICE1 was reversely altered in OXs and KO, pointing to a potential cause for the observed changes in stomatal phenotypes. Stomatal and transcriptional changes were not observed in the Arabidopsis rfs double mutant. Taken together, RFS is a novel regulator of SD and SS and is a promising candidate for genetic engineering of climate-resilient crops.

The reference genome and abiotic stress responses of the model perennial grass Brachypodium sylvaticum.

Perennial grasses are important forage crops and emerging biomass crops and have the potential to be more sustainable grain crops. However, most perennial grass crops are difficult experimental subjects due to their large size, difficult genetics, and/or their recalcitrance to transformation. Thus, a tractable model perennial grass could be used to rapidly make discoveries that can be translated to perennial grass crops. Brachypodium sylvaticum has the potential to serve as such a model because of its small size, rapid generation time, simple genetics, and transformability. Here, we provide a high-quality genome assembly and annotation for B. sylvaticum, an essential resource for a modern model system. In addition, we conducted transcriptomic studies under 4 abiotic stresses (water, heat, salt, and freezing). Our results indicate that crowns are more responsive to freezing than leaves which may help them overwinter. We observed extensive transcriptional responses with varying temporal dynamics to all abiotic stresses, including classic heat-responsive genes. These results can be used to form testable hypotheses about how perennial grasses respond to these stresses. Taken together, these results will allow B. sylvaticum to serve as a truly tractable perennial model system.

Characterization of diterpene synthase genes in Brachypodium distachyon, a monocotyledonous model plant, provides evolutionary insight into their multiple homologs in cereals.

Gibberellins are diterpenoid phytohormones that regulate plant growth, and are biosynthesized from a diterpene intermediate, ent-kaurene, which is produced from geranylgeranyl diphosphate via ent-copalyl diphosphate (ent-CDP). The successive 2 cyclization reactions are catalyzed by 2 distinct diterpene synthases, ent-CDP synthase (ent-CPS) and ent-kaurene synthase (KS). Various diterpene synthase genes involved in specialized metabolism were likely created through duplication and neofunctionalization of gibberellin-biosynthetic ent-CPS and KS genes in crops. Brachypodium distachyon is a monocotyledonous species that is a model plant in grasses. We herein found 1 ent-CPS gene homolog BdCPS and 4 tandemly arrayed KS-like genes BdKS1, KSL2, KSL3, and KSL4 in the B. distachyon genome, a simpler collection of paralogs than in crops. Phylogenetic and biochemical analyses showed that BdCPS and BdKS1 are responsible for gibberellin biosynthesis. BdKSL2 and BdKSL3 are suggested to be involved in specialized diterpenoid metabolism. Moreover, we restored KS activity of BdKSL2 through amino acid substitution.

PHYTOCHROME C regulation of photoperiodic flowering via PHOTOPERIOD1 is mediated by EARLY FLOWERING 3 in Brachypodium distachyon.

Daylength sensing in many plants is critical for coordinating the timing of flowering with the appropriate season. Temperate climate-adapted grasses such as Brachypodium distachyon flower during the spring when days are becoming longer. The photoreceptor PHYTOCHROME C is essential for long-day (LD) flowering in B. distachyon. PHYC is required for the LD activation of a suite of genes in the photoperiod pathway including PHOTOPERIOD1 (PPD1) that, in turn, result in the activation of FLOWERING LOCUS T (FT1)/FLORIGEN, which causes flowering. Thus, B. distachyon phyC mutants are extremely delayed in flowering. Here we show that PHYC-mediated activation of PPD1 occurs via EARLY FLOWERING 3 (ELF3), a component of the evening complex in the circadian clock. The extreme delay of flowering of the phyC mutant disappears when combined with an elf3 loss-of-function mutation. Moreover, the dampened PPD1 expression in phyC mutant plants is elevated in phyC/elf3 mutant plants consistent with the rapid flowering of the double mutant. We show that loss of PPD1 function also results in reduced FT1 expression and extremely delayed flowering consistent with results from wheat and barley. Additionally, elf3 mutant plants have elevated expression levels of PPD1, and we show that overexpression of ELF3 results in delayed flowering associated with a reduction of PPD1 and FT1 expression, indicating that ELF3 represses PPD1 transcription consistent with previous studies showing that ELF3 binds to the PPD1 promoter. Indeed, PPD1 is the main target of ELF3-mediated flowering as elf3/ppd1 double mutant plants are delayed flowering. Our results indicate that ELF3 operates downstream from PHYC and acts as a repressor of PPD1 in the photoperiod flowering pathway of B. distachyon.

High-resolution Hi-C maps highlight multiscale chromatin architecture reorganization during cold stress in Brachypodium distachyon.

BACKGROUND: The adaptation of plants to cold stress involves changes in gene expression profiles that are associated with epigenetic regulation. Although the three-dimensional (3D) genome architecture is considered an important epigenetic regulator, the role of 3D genome organization in the cold stress response remains unclear. RESULTS: In this study, we developed high-resolution 3D genomic maps using control and cold-treated leaf tissue of the model plant Brachypodium distachyon using Hi-C to determine how cold stress affects the 3D genome architecture. We generated ~ 1.5 kb resolution chromatin interaction maps and showed that cold stress disrupts different levels of chromosome organization, including A/B compartment transition, a reduction in chromatin compartmentalization and the size of topologically associating domains (TADs), and loss of long-range chromatin loops. Integrating RNA-seq information, we identified cold-response genes and revealed that transcription was largely unaffected by the A/B compartment transition. The cold-response genes were predominantly localized in compartment A. In contrast, transcriptional changes are required for TAD reorganization. We demonstrated that dynamic TAD events were associated with H3K27me3 and H3K27ac state alterations. Moreover, a loss of chromatin looping, rather than a gain of looping, coincides with alterations in gene expression, indicating that chromatin loop disruption may play a more important role than loop formation in the cold-stress response. CONCLUSIONS: Our study highlights the multiscale 3D genome reprogramming that occurs during cold stress and expands our knowledge of the mechanisms underlying transcriptional regulation in response to cold stress in plants.

Submergence Stress Alters the Expression of Clock Genes and Configures New Zeniths and Expression of Outputs in Brachypodium distachyon.

Plant networks of oscillating genes coordinate internal processes with external cues, contributing to increased fitness. We hypothesized that the response to submergence stress may dynamically change during different times of the day. In this work, we determined the transcriptome (RNA sequencing) of the model monocotyledonous plant, Brachypodium distachyon, during a day of submergence stress, low light, and normal growth. Two ecotypes of differential tolerance, Bd21 (sensitive) and Bd21-3 (tolerant), were included. We submerged 15-day-old plants under a long-day diurnal cycle (16 h light/8 h dark) and collected samples after 8 h of submergence at ZT0 (dawn), ZT8 (midday), ZT16 (dusk), ZT20 (midnight), and ZT24 (dawn). Rhythmic processes were enriched both with up- and down-regulated genes, and clustering highlighted that the morning and daytime oscillator components (PRRs) show peak expression in the night, and a decrease in the amplitude of the clock genes (GI, LHY, RVE) was observed. Outputs included photosynthesis-related genes losing their known rhythmic expression. Up-regulated genes included oscillating suppressors of growth, hormone-related genes with new late zeniths (e.g., JAZ1, ZEP), and mitochondrial and carbohydrate signaling genes with shifted zeniths. The results highlighted genes up-regulated in the tolerant ecotype such as METALLOTHIONEIN3 and ATPase INHIBITOR FACTOR. Finally, we show by luciferase assays that Arabidopsis thaliana clock genes are also altered by submergence changing their amplitude and phase. This study can guide the research of chronocultural strategies and diurnal-associated tolerance mechanisms.

Analysis of Brachypodium distachyon UVR8 reveals conservation in UV-B receptors.

The Ultraviolet Resistance Locus 8 (UVR8) in plants recognizes ultraviolet-B (UV-B) light and plays a crucial role in regulating plant growth through a series of signal transduction events. However, the UVR8 in monocotyledon crops has not yet been systematically analysed. We identified BdUVR8 (BRADI_3g45740) from the genome of Brachypodium distachyon, a relative of wheat, by analysing the phylogenetic tree, the gene expression pattern, detecting accumulation of UV-B response metabolites, and checking for phenotype recovery. The BdUVR8 protein sequence is similar to the known UVR8 of other species. The phylogenetic tree of UVR8 shows clear divergence between dicotyledons and monocotyledons. Expression analysis revealed that UV-B downregulates BdUVR8 by 70% and upregulates the chalcone synthase (BdCHS) gene 3.4-fold in B. distachyon. The pCAMBIA1300::BdUVR8-mCherry construct introduced into Arabidopsis uvr8 mutants showed that the BdUVR8 protein is localized in the cytoplasm and translocates into the nucleus in response to UV-B irradiation. The introduction of BdUVR8 into uvr8 rescued hypocotyl elongation caused by UV-B and restored expression of HY5, Chalcone synthase, and Flavanone 3-hydroxylase, as well as accumulation of total flavonoids. Together, our results show that BdUVR8 is a photoreceptor that perceives UV-B in B. distachyon.

Regulation of hair cell and stomatal size by a hair cell-specific peroxidase in the grass Brachypodium distachyon.

The leaf epidermis is the outermost cell layer forming the interface between plants and the atmosphere that must both provide a robust barrier against (a)biotic stressors and facilitate carbon dioxide uptake and leaf transpiration.1 To achieve these opposing requirements, the plant epidermis developed a wide range of specialized cell types such as stomata and hair cells. Although factors forming these individual cell types are known,2,3,4,5 it is poorly understood how their number and size are coordinated. Here, we identified a role for BdPRX76/BdPOX, a class III peroxidase, in regulating hair cell and stomatal size in the model grass Brachypodium distachyon. In bdpox mutants, prickle hair cells were smaller and stomata were longer. Because stomatal density remained unchanged, the negative correlation between stomatal size and density was disrupted in bdpox and resulted in higher stomatal conductance and lower intrinsic water-use efficiency. BdPOX was exclusively expressed in hair cells, suggesting that BdPOX cell-autonomously promotes hair cell size and indirectly restricts stomatal length. Cell-wall autofluorescence and lignin stainings indicated a role for BdPOX in the lignification or crosslinking of related phenolic compounds at the hair cell base. Ectopic expression of BdPOX in the stomatal lineage increased phenolic autofluorescence in guard cell (GC) walls and restricted stomatal elongation in bdpox. Together, we highlight a developmental interplay between hair cells and stomata that optimizes epidermal functionality. We propose that cell-type-specific changes disrupt this interplay and lead to compensatory developmental defects in other epidermal cell types.

Complexity and regulation of age-dependent alternative splicing in Brachypodium distachyon.

Alternative splicing (AS) is a gene regulatory mechanism that generates multiple transcripts of the same gene precursor by the spliceosome complex, promoting messenger RNA complexity, and proteome diversity. Although AS is extensively studied in response to environmental stresses, whether it mediates age-dependent development and how it is adjusted by growth transitions are largely unknown. Here, we comprehensively explored the AS landscape at different developmental stages in the grass model plant Brachypodium (Brachypodium distachyon). We identified abundant coding genes and noncoding transcripts subject to dynamic AS regulation during juvenile, adult, and reproductive transitions. Moreover, we revealed that SC35-LIKE SPLICING FACTOR 33 (SCL33), a serine/arginine-rich splicing factor in spliceosomes, plays a redundant and antagonistic role with its putative paralog, SCL33L, in regulating intron assembly across distinct developmental stages. In addition, we determined global AS variations in microRNA156 (miR156)-overproducing plants, in which growth transitions are delayed, and found that SPLs were regulated by miR156 in intron retention alteration in addition to mRNA clearance and translation inhibition manners. Finally, we demonstrated a complex regulatory process of age-dependent AS events in B. distachyon that was coincidently or separately regulated by miR156 and SCL33/SCL33L. These results illustrate a substantial machinery of AS that mediates phase transitions in plants.

Mixed-Linkage Glucan Is the Main Carbohydrate Source and Starch Is an Alternative Source during Brachypodium Grain Germination.

Seeds of the model grass Brachypodium distachyon are unusual because they contain very little starch and high levels of mixed-linkage glucan (MLG) accumulated in thick cell walls. It was suggested that MLG might supplement starch as a storage carbohydrate and may be mobilised during germination. In this work, we observed massive degradation of MLG during germination in both endosperm and nucellar epidermis. The enzymes responsible for the MLG degradation were identified in germinated grains and characterized using heterologous expression. By using mutants targeting MLG biosynthesis genes, we showed that the expression level of genes coding for MLG and starch-degrading enzymes was modified in the germinated grains of knocked-out cslf6 mutants depleted in MLG but with higher starch content. Our results suggest a substrate-dependent regulation of the storage sugars during germination. These overall results demonstrated the function of MLG as the main carbohydrate source during germination of Brachypodium grain. More astonishingly, cslf6 Brachypodium mutants are able to adapt their metabolism to the lack of MLG by modifying the energy source for germination and the expression of genes dedicated for its use.